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. Abstract 
IT} . _ 

^-j- ' We report on a search for D° — D° mixing made by a study of 

^— l ■ the 'wrong-sign' process D° — > K + ir~ . The data come from an in- 

tegrated luminosity of e + e~ collisions at y/s m 10 GeV consisting 

> 1 of 9.0 ft)" 1 , recorded with the CLEO-II.V detector. We measure 

the time-integrated rate of the 'wrong-sign' process 
. relative to that of the Cabibbo- favored process D° — > K + ir~ , to be 

Q\ ; R ws = (0.34 ± 0.07 ± 0.06)%. We study that rate as a function of 

O ■ the decay time of the D°, to distinguish the rate of direct doubly- 

. Cabibbo-suppressed decay from D° — D° mixing. The amplitudes 

that describe D° — D° mixing, x' and y', are consistent with zero. 
• The one-dimensional limits, at the 95% C.L., that we determine are 

jL^: (l/2)x /2 < 0.05%, and -5.9% < y' < 0.3%. All results are prelimi- 

D : nary. 

^ : 

> : Studies of the evolution of a K° or B® into the respective anti-particle, 

a K° or Bj, have guided the form and content of the Standard Model, 
and permitted useful estimates of the masses of the charm and top quark 
masses, prior to direct observation of those quarks. In this paper, we present 
the results of a search for the evolution of the D° into the D°, where the 
principal motivation is to glimpse new physics outside the Standard Model, 
prior to direct observation of that physics at the high energy frontier. 

A D° can evolve into a D° through on-shell intermediate states, such as 
K + K~ with rriK+K- =^d°, or through off-shell intermediate states, such 
those that might be present due to 'new physics'. We denote the amplitude 
through the former (latter) states by — iy (x), in units of T D o/2 0, and we 
neglect all types of CP violation. 

For comparison, in the K° — K° system the analogous y and \x\ are 
both near unity 0, |3|. The prediction that an initial K° will decay with 
two lifetimes, as described by the non-zero y, was famously made by Gell- 
Mann and Pais 0. In the B^ — B Q d system theory firmly predicts that y 
is negligible, and experiments have not yet sought y out. For the D° — D° 
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system, y is likely to receive significant contributions from the Standard 
Model ||. It may even be that y dominates the D° — > D° amplitude. 

It has been x, for both the K° — K° and B® — B° d systems, that has 
provided information about the charm and top quark masses. A report of 
the first measurement in 1961 of \x\ for the K°—K° system noted 'We cannot 
compare our experimental value for \x\ with any theoretical calculation' ||. 
By 1974, theory had caught up, and exploited \x\ for the K° — K° system 
to predict the charm quark mass [0], just before that quark was discovered. 
The tiny, CP violating Im(x) for the K° is sensitive to the value of the top 
quark mass. The large value for the B% of \x\ (now 0.73±0.03 0) indicated 
that the top quark is very massive ||, in distinction to contemporaneous 
data from the high energy frontier ||. 

Many predictions for x in the D° — > D° amplitude have been made [[10 



The Standard Model contributions are suppressed down to at least \x\ ~ 
tan 2 Q c 0.05 because D° decay is Cabibbo-favored; the GIM |TTJ cancel- 



lation could suppress \x\ down to 10 —0.01. Many non-Standard Models, 
particularly those that address patterns of quark flavor, predict \x\ ~ 0.01 
or greater. Contributions to x at this level can result from the presence 



of new particles with masses as high as 100 TeV [I2|. Decisive signatures 
of such particles might include \y\ <C \x\, or CP-violating interference be- 
tween a substantial imaginary component of x and either y, or a direct 
decay amplitude. In order to accurately assess the origin of a D°—D° mix- 
ing signal, the effects described by y must be distinguished from those that 
are described by x. 

We report here on a study of the process D° — > K + tt~ . We use the 
charge of the 'slow' pion, 7r+, from the decay D* + — ► -D°7r+ to deduce pro- 
duction of the D , and then we seek the rare 'wrong-sign' K + tt~ final state 
(WS), in addition to the more frequent 'right-sign' final state, K~tc + (RS). 
The wrong-sign process, D° — > K + tt~ , can proceed either through direct 
doubly- Cabibbo-suppressed decay (DCSD), or through mixing followed by 
the Cabibbo-favored decay (CFD), D° — ► D° — > K + n~. Both processes 
contribute to the time integrated 'wrong-sign' rate, i? ws : 

T(D°^K+it-) 



To disentangle the two processes that could contribute to D° — > K + n~ , 
we study the distribution of wrong-sign final states as a function of the 
proper decay time, t, of the D°. The proper decay time is in units of the 
mean D° lifetime, r^o = 415 ± 4fs The differential 'wrong-sign' rate, 
relative to T(D° ^K + 7r~), is r ws (t) 0: 



s (t) = [R D + y^y't + I(x /2 +2 / 2 )t 2 ]e-*, (1) 
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where the modified mixing amplitudes x' and y' in Eqn. |I] are given by: 



y' = y cos 5— x sin 5 
x' = x cos 5+y sin 5 ' 

and 5 is a possible strong phase between DCSD and CFD amplitudes: 

-jR D e~ 



,-it - (K + ir-\T\D°) 



(K+7T-\T\D°) 

The symbol Re> represents the DCSD rate, relative to the CFD rate. There 
are plausible theoretical arguments that S is small Jl5, [l(J. If both D° and 



D° exclusively populate the I = 1/2 amplitude of K + n , then 5 would 
be zero. For the CFD, D° -> K+n-, the 1 = 3/2 amplitude is indeed 
disfavored, with 1^43/2/^1/2! = 0.27 ±0.02 |I7|] . A non-zero 6 could develop 
if D° — > K + 7i~ populates the 1 = 3/2 amplitude differently than the CFD 
does. Crudely, one might guess \5\ < ~ | ^3/2/^1/2 1 ~ 15°. The size of S could 
be settled by measurements of the DCSD contributions to D + — > K + tt°, 
D + ^K^ti + , and D° — >K^it , which are now feasible with our data set. 

For decays to wrong-sign final states other than K + ir~ , such as K + tt~7t°, 
or if + 7r~7r + 7r~, there will be distinct strong phases. Moreover, the broad 
resonances that mediate those multibody hadronic decays, such as the K*, 
p, etc., modulate those phases as a function of position on the Dalitz plot. 
Thus, multibody hadronic wrong-sign decays might afford an opportunity 
to distinguish x and y from x' and y'. 

An important aspect of the expression Eqn. |I| for the two-body decay 
D° — > K + n~ is that the dependence on y' and x' is distinguishable due to 
the interference with the direct decay amplitude, which induces the term 
linear in t. Such behavior is complementary to the differential decay rate 
to CP eigenstates such as D° — > K + K~, which is sensitive to y alone, or 
that of D° ^K + £~V e , which is sensitive to r mix = (x 2 +y 2 )/2 = (x' 2 +y' 2 )/2 
alone. 

Our data was accumulated between Jan. 1996 and Feb. 1999 from an in- 
tegrated luminosity of e + e~ collisions at y's ~ 10 GeV consisting of 9.0 fb^ 1 , 
provided by the Cornell Electron Storage Rings (CESR). The data were 
taken with CLEO-II multipurpose detector | 18|| , upgraded in 1995 when a 



silicon vertex detector (SVX) was installed [0 and the drift chamber gas 
was changed from argon-ethane to helium-propane. The upgraded config- 
uration is named CLEO-II. V, where the 'V is short for Vertex. 

We reconstruct candidates for the decay sequences D* + — > vr+Z)°, fol- 
lowed by either D° — >K + tt~ (wrong-sign) or 7r + (right-sign). The 
sign of the slow charged pion tags the charm state at t = as either D° or 
D° . The broad features of the reconstruction are similar to those employed 
in the recent CLEO-II. V measurement of the D meson lifetimes [EUI. The 
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following points are important in understanding how we have improved our 
sensitivity to D° — > K + n~ , relative to earlier CLEO work pTl EEJ: 



The SVX allows substantially more precise measurement of charged 
particle trajectories in the dimension parallel to the colliding beam 
axis. When combined with improvements in track-fitting, the CLEO- 
II. V resolution for Q = M^ s ktt ~ M, where M is the mass of the 
K + 7t~ system, is ctq = 190 ± 6 KeV, compared with the earlier value 
of a Q = 750 KeV ^Tj. 



2. The use of helium-propane, in addition to improvements in track- 
fitting, have reduced the CLEO-II.V resolution for M to <jm = 6.4 ± 



0.1 MeV, compared with the earlier value of <jm = 11 MeV [21 



3. Improved mass resolution, as well as rejections based on the mo- 
mentum asymmetry of the two charged tracks, allow clean separa- 
tion of the signal D° — > K + n~ from D° — » K + K~ , D° — ► , and 
D° — ► K~7i + , and from multibody decays of the D°, at a cost of about 
35% of the signal acceptance. We use the modest ix + / K + separation 
provided by measurement of dE/dx in the CLEO-II.V drift chamber 
primarily for systematic studies. Addition of a new device with per- 
fect tt + /K + identification and acceptance might enable the recovery 
of the 35% acceptance loss. 

Multiple scattering on the field wires, which constitute 70% of the ma- 
terial, as measured in radiation lengths, in the CLEO-II.V drift chamber, 
appears to dominate the current ctq and <jm- Should the track-fitting be 
altered to treat scattering on the field wires as discrete and localized, both 
ctq and om might improve by as much as a factor of 2. 

Our signal for the wrong-sign process D° — > K + ir~ is shown in Fig. [I]. 
We determine the background levels by performing a fit to the plane of 
0<Q<20 MeV versus 1.76 <M< 1.97 GeV, which has an area about 150 
times larger than our signal region. Event samples generated by the Monte- 
Carlo method and fully simulated in our detector, corresponding to 90 fb _1 
of integrated luminosity, are used to estimate the background shapes in the 
Q—M plane. The shapes are allowed to float in a fit to the data; the results 
of the fit are displayed in Fig |1[ The excess of events in the signal region 
is prominent. 

We describe the signal shape with the right-sign data that is within 7 a 
of the nominal CFD value in the Q—M plane. The results of the fit to the 
wrong-sign data are summarized in Table I. 

No acceptance corrections are needed to compute R ws = (0.34 ± 0.07)% 
from Table I. The dominant systematic errors all stem from the potentially 
inaccurate modeling of the initial and acceptance-corrected shapes of the 
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Figure 1: Signal for the wrong-sign process D° — > K + n~ . For the top plot, 
M is within 14 MeV of the nominal CFD value, and for the bottom plot, Q 
is within 500 KeV of the nominal CFD value. The data are the full circles 
with error bars, the fit to the signal is cross-hatched, and the fits to the 
backgrounds are singly hatched. The results of the fit are summarized in 
Table I. 
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background contributions in the Q-M plane. We assess these systematic 
errors by substantial variation of the fit regions, dE/dx criteria, and kine- 
matic criteria; the total systematic error we assess is 0.06%. 

Our complete result for R ws is summarized in Table II. As we describe 
later, our data are consistent with an absence of D° — D° mixing, so our 
best estimate for the relative DCSD rate, Rjj, is that it equals i? ws . 

There are two directly comparable measurements of R ws : one is from 
CLEO-II HI], R ws = (0.77 ± 0.25 ± 0.25)% which used a data set indepen- 
dent of that used here; the second is from Aleph |23|, R ws = (1.84 ± 0.59 ± 



0.34)%; comparison of our result and these are marginally consistent with 



X 



6.0 for 2 DoF, for a C.L. of 5.0%. 



We have split our sample into candidates for D° — > K + n~ and D° — > 
K~7i + . There is no evidence for a CP-violating time-integrated asymmetry. 
From Table I, it is straightforward to estimate the la statistical error on 
the CP violating time-integrated asymmetry as V107/54.8 ~ 20%. 

TABLE I. Event yields in a signal region of 2.4 a centered on the nominal 
Q and M values. The total number of candidates is 107. The bottom row 
describes the normalization sample. 



Component 


# Events 


D°^K + 7i- (WS Signal) 


54.8 ± 10.8 


random n ± + D°/D° 


24.3 ± 1.8 


cc 


12.3 ±0.8 


uds 


8.6 ±0.4 


D° — > Pseudoscalar- Vector 


7.0 ±0.4 



D° -> K + 7v- (RS Normalization) 16126 ± 126 



TABLE II. Result for R ws . For the branching ratio B(D° — > K + tt~) we take 
the absolute branching ratio B(D° -^K + tt~) = (3.85±0.09)%. The third error 
results from the uncertainty in this absolute branching ratio. As discussed in 
the text, our best estimate is that Rd = Rws- 

Quantity Result 

R^ (0.34 ±0.07 ±0.06)% 

R ws /tan 4 9c (1.28 ± 0.25 ± 0.21) 

B(D ^K + 7T-) (1.31 ± 0.26 ± 0.22 ± 0.03) x 10~ 4 

Given the absence of a significant time-integrated CP asymmetry, we 
undertake a study of the decay time dependence wrong-sign rate based 
upon Eqn. |I[ We reconstruct the proper D° decay time primarily from the 
vertical displacement of the K + tt~ vertex from the e + e~ collision 'ribbon,' 
which is infinitesimal in its vertical extent. We require a well-reconstructed 
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vertex in 3- dimensions, which causes a loss of about 12% of the candi- 
dates described in Table I. Our resolution, in units of the mean D° life, 
is about 1/2. Study of the plentiful right-sign sample allows us to fix our 
detailed resolution function, and shows that biases in the reconstruction of 
the proper decay time contribute negligibly to the wrong-sign results. 




I I I I I 

-2.5 2.5 5 7.5 

D° Lifetimes 



Figure 2: Distribution in t for the D — > K + -k~ candidates. The data are 
shown as the full circles with error bars. All other information comes from 
the fit to the data, The smooth curves show the various specific contri- 
butions as labeled. The cross hatched region is the net contribution from 
D° — > K + n~ , after incorporation of the (destructive) interference and mix- 
ing. 

The distribution of proper decay times, t, for wrong-sign candidates that 
are within 2.4 a of the nominal CFD value in the Q — M plane is shown in 
Fig. [2| Maximum-likelihood fits are made to those data. The backgrounds 
are described by levels and shapes deduced from the fit to the Q—M plane, 
and from study of the simulated data sample. The wrong-sign signal is 
described by Eqn. p], folded with our resolution function. 



7 



For the benchmark fit to the wrong-sign data, x' and y' are constrained 
to be zero. This fit has a confidence level of 84%, indicating a good fit. 

The mixing amplitudes x' and y' are then allowed to freely vary, and the 
best fit values are shown in Fig. || and in Table III. The fit improves slightly 
when x' and y' are allowed to float to the values that maximize the likeli- 
hood. However, the small value of the likelihood change a/~2A ln£ = 1.8 o 
(including systematic errors) does not permit us to eliminate the possibility 
that the improvement is due to a statistical fluctuation. 

Therefore, our principal results concerning mixing are the one-dimensional 
intervals, which correspond to a 95% confidence level, that are given in the 
second column of Table III. 

Additionally, we evaluate a contour in the two-dimensional plane of y' 
versus x' which, at 95% confidence level, contains the true value of x' and 
y'. To do so, we determine the contour around our best fit values where 
the — ln£ has increased by 3.0 units. All other fit variables, including the 
DCSD rate and background contributions, are allowed to float to distinct, 
best fit values at each point on the contour. The interior of the contour is 
shown, as the small, dark, cross-hatched region near the origin of Fig. |3|. On 
the axes of x' and y', this contour falls slightly outside the one- dimensional 
intervals listed in Table III, as expected. 

TABLE III. Results of the fit, where x' and y' are free to float, to the distri- 
bution of D° — > K + ti~ candidates in t. The data and the fit components are 
shown in Fig. ^|. 



Parameter 


Best Fit 


95% C.L. 


Rd 

y' 


(0.501^2 ± 0.08)% 
(-2.71U ±0.2)% 


0.22% <R D < 0.77% 
-5.9% < y' < 0.3% 


x' 

(l/2)x' 2 


(0± 1.6 ±0.2)% 


\x'\ < 3.2% 
< 0.05% 



We have evaluated the allowed regions of other experiments p4| , Pq , |26 
p7| at 95% C.L., and shown those regions in Fig. 
All results described here are preliminary. 

If we assume that 5 is small, which is plausible [JO], [RJ then x' ~x and 
we can indicate the impact of our work in limiting predictions of D° — D° 
mixing from extensions to the Standard Model. Eighteen of the predictions 



recently tabulated |10| have some inconsistency with our limit. Among 
those predictions, some authors have made common assumptions, however. 

Because our data is consistent with an absence of D° — D° mixing, our 
best information on the DCSD rate, Rd, is that it equals R ws , as summa- 
rized in Table II. 

We will soon complete studies where we allow various types of CP 
violation to modify Eqn. [1]. Also, the D° decay modes K + K~, 7r + 7r~, 
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Figure 3: Limits in the y' vs. x' plane. Our experiment limits, at 95% 
C.L., the true values of x' and y' to occupy the cross-hatched region near the 
origin. Also shown are the similar zones from other recent experiments. We 
assume 5 — to place the recent work of E791 that utilized D° — > K + K~; 
a non-zero 5 would rotate the E791 confidence region clockwise about the 
origin by an angle of 5. 
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K$<f), Ksir + iT~, K + 7t~7t°, K + 7t~7t + 7t~ , and K + l~Vi are under active in- 
vestigation, using the CLEO-II.V data. The (Ksir + ) resonance bands in 
D° — > Ks7t + tt~ permit investigation of D°^D°, with sensitivity heightened 
by coherent interference that is modulated by resonant phases. 

The entire CLEO-II.V data set, suitably exploited, could be used to 
observe D° — D° mixing if either \x\ or \y\ exceed approximately 1%. 
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